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Abstract 



In this paper we study the computational complexity of computing the noncommutative 
determinant. We first consider the arithmetic circuit complexity of computing the noncom- 
mutative determinant polynomial. Then, more generally, we also examine the complexity of 
computing the determinant (as a function) over noncommutative domains. Our hardness re- 
sults are summarized below: 

• We show that if the noncommutative determinant polynomial has small noncommuta- 
tive arithmetic circuits then so does the noncommutative permanent. Consequently, the 
commutative permanent polynomial has small commutative arithmetic circuits. 

• For any field F we show that computing the n x n permanent over F is polynomial-time 
reducible to computing the 2n x 2?! (noncommutative) determinant whose entries are 
0(71^) X 0{n^) matrices over the field F. 

• We also derive as a consequence that computing the n x n permanent over nonnegative 
rationals is polynomial-time reducible to computing the noncommutative determinant over 
Clifford algebras of n'-"^^^ dimension. 

Our techniques are elementary and use primarily the notion of the Hadamard Product of non- 
commutative polynomials. 



1 Introduction 



In his seminal paper |N91) Nisan first systematically studied the problem of proving lower bounds 
for noncommutative computation. The focus of his study was noncommutative arithmetic circuits, 
noncommutative arithmetic formulas and noncommutative algebraic branching programs. In his 
central result based on a rank argument, Nisan shows that the noncommutative permanent or deter- 
minant polynomials in the ring ¥{xii, • • • , require exponential size noncommutative algebraic 
branching programs. 

Nisan's results are over the free noncommutative ring ¥{X). Chien and Sinclair, in |CS04j . explore 
the same question over other noncommutative algebras. They refine Nisan's rank argument to show 
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exponential size lower bounds for formulas computing the permanent or determinant over specific 
noncommutative algebras, like the algebra of 2 x 2 matrices over F, the quaternion algebra, and a 
host of other examples. 

However, the question of whether there is a small noncommutative circuit for the determinant 
or permanent remains unanswered. (Indeed, no explicit lower bounds are known for the general 
noncommutative circuit model.) Since the existence of small noncommutative arithmetic circuits 
for the permanent would imply the existence of small commutative arithmetic circuits for the per- 
manent, we have a good reason to believe that the permanent does not have small noncommutative 
arithmetic circuits. However, as far as we know, no such argument has been given for the case of 
the noncommutative determinant. Indeed, since Nisan ^N9lj has also shown an exponential sepa- 
ration between the power of noncommutative formulas and circuits, it may very well be that the 
noncommutative determinant has polynomial-sized arithmetic circuits. 

Another motivation for studying the computational difficulty of computing the noncommutative 
determinant (as a function) is an approach to designing randomized approximation algorithms 
for the — 1 permanent by designing good unbiased estimators based on the determinant. This 
approach has a long history starting with [GGSlj |KKL+93] . Of specific interest are the works 
of Barvinok [B]; Chien, Rasmussen, and Sinclair [CRS03] : and more recently that of Moore and 
Russell |MR09) . Barvinok [B] defines a variant of the noncommutative determinant called the 
symmetrized determinant and shows that given inputs from a constant dimensional matrix alge- 
bra, the symmetrized determinant over these inputs can be evaluated in polynomial time. He uses 
these to define a series of algorithms that he conjectures might yield progressively better random- 
ized approximation algorithms for the (commutative) permanent. Chien, Rasmussen, and Sinclair 
|CRS03] show that efficient algorithms to compute the determinant over Clifford algebras of poly- 
nomial dimension would yield efficient approximation algorithms for the permanent. Moore and 
Russell |MR09| provide evidence that Barvinok's approach might not work, but their results also 
imply that computing the symmetrized or standard noncommutative determinant over polynomial 
dimensional matrix algebras would give a good estimator for the permanent. 

Our results 

1. We provide evidence that the noncommutative determinant is hard. We show that if the 
noncommutative determinant^ can be computed by a small noncommutative arithmetic cir- 
cuit, then so can the noncommutative permanent and therefore, the commutative permanent 
has small commutative arithmetic circuits. This is in marked contrast to the commutative 
case, where the determinant is known to be computable by polynomial sized circuits, but the 
permanent is not known (or expected) to have subexponential sized arithmetic circuits. 

2. We show that computing the noncommutative determinant over matrix algebras of polynomial 
dimension is as hard as computing the commutative permanent. We also derive as a conse- 
quence that computing the n x n permanent over nonnegative rationals is polynomial-time 
reducible to computing the noncommutative determinant over Clifford algebras of poly(n) 
dimension. 

^We haven't defined this polynomial formally yet and there are, in fact, many ways of doing it. See Section (2] 
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This points to the intractabiUty of carrying over Barvinok's approach for large dimension, 
and also to the possibility that the approach of Chien, Rasmussen, and Sinclair might be 
computationally infeasible. 

We stress that our result here is potentially more useful than a noncommutative circuit lower 
bound for the determinant, from an algorithmic point of view. For, an arithmetic circuit 
lower bound result would not rule out the possibility of a polynomial-time algorithm for the 
noncommutative determinant over even polynomial dimension matrix algebras. For example, 
Barvinok's algorithm [BJ computes the symmetrized determinant over constant dimensional 
matrix algebras, whereas any algebraic branching program that computes the symmetrized 
determinant over constant dimensional matrix algebras must be of exponential size |CS04| . 

2 Preliminaries 

For any set of variables X, let ¥{X) denote the ring of noncommuting polynomials over X. Let 
Ai{X) denote the set of noncommutative monomials over X; given d S N, let Aiai^) denote the 
monomials over X of degree exactly d. For / € F{X) and m G Ai{X), we will denote by f{m) the 
coefficient of the monomial m in /. 

For any ring R, we use Mn{R) to denote the ring of n x n matrices with entries from R. 

Fix X = {xi, X2, . . . , Xm} and Y = {yi, y2, ■ ■ ■ , Un}, two disjoint sets of variables. Given / G F(X), 
matrices Ai € Mk{¥{Y)) for 1 < i < m, and io,io ^ [k], we use f{Ai,A2, ■ ■ ■ , Am){io, jo) to denote 
the (zo, jo)th entry of the matrix f{Ai,A2, . . . , Am) € Mk{¥{Y)). 

2.1 Noncommutative determinants and permanents 

Given X = {xij | 1 < i, j < ?^} for n G N, we define the n x n noncommutative determinant and 
permanent polynomials over the set of variables X. By fixing the order of multiplication in each 
monomial of the commutative determinant /permanent polynomials in different ways, one can obtain 
many different reasonable ways of defining the nx n noncommutative determinant and permanent, 
and indeed many of these definitions have been studied (see |A96] . which surveys various flavours of 
the noncommutative determinant). The most straightforward definitions are those of the Cayley de- 
terminant and Cayley permanent - we will denote these by Cdetri(X) and Cperm„(X) respectively 
- which use the row order of multiplication. That is, 

Cdet„(X) = sgn(f^) 
Cperm„(X) = ^ 

2^1, (t(1) ■ ^2,(7(2) ■ ■ ■ Xn,cr{n)- 

We also define the Moore determinant and Moore permanent - denoted Mdet„(X) and Mperm^(X) 
respectively - by ordering the variables in each monomial using the cyclic order of the 
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corresponding permutation. Given a G /S^, we write it as a product of disjoint cycles 
{nl^ ■ ■ ■ n'[i^)(n2i ■ ■ ■ n2i^) ■ ■ ■ {n'^^ ■ ■ ■ n"i^) such that Mi € [r] and j G \lr\ \ {1}, we have n'^^ < n"- and 
71^]^ > > • • • > n"i. The Moore determinant and permanent are defined as 



Mdet„(X) = ^ sgn((T) x^j^.n^a ' • • ^"y,,"?! ' ' ' ^Ki,K2 ' ' ' ^< 

<J€Sn 

Mperm„(X) = ^ a;„j^,„j^ • • • • • • x^^^^n^^ • • • Xn-^ 



<^n", ,n", ■ 

crGSn 



In the setting of a field F of characteristic 0, Alexander Barvinok, in [Bj, has studied another 
variant of the noncommutative determinant called the symmetrized determinant, which is denoted 
sdetn(X). It is defined as follows: 

sdet„(X) = X] Sgn(o-)sgn(T) X^(i)_^(i)X^(2),^(2) • • • X^(n),a(n)- 



Barvinok shows that, for any fixed dimensional associative algebra A over F of characteristic zero, 
there is a polynomial-time algorithm which, on input an n x n matrix A with entries from A, 
computes sdetn(A). It is not known whether such algorithms exist for the Cayley or Moore deter- 
minants. 



2.2 Models for noncommutative arithmetic computation 

A noncommutative arithmetic circuit C over a field F is defined as follows: C is a directed acyclic 
graph and every leaf of the graph is labelled with either an input variable from the set of variables 
X or an element from F. Every internal node is labelled by either (-I-) or (x) - meaning that it is 
either an addition or multiplication gate respectively - and has fanin two. Since we are working 
over noncommutative domains, we will assume that each multiplication gate has a designated left 
child and a designated right child. Each gate of the circuit computes a polynomial in ¥{X) in the 
natural way: the polynomials computed at the leaves are the polynomials labelling the leaves; the 
polynomial computed at an internal node labelled by -|- (resp. x) is the sum (resp. product in 
left-to-right order) of the polynomials computed at its children. The polynomial computed by C is 
the polynomial computed at a designated output node of the circuit. 

We also recall the definition of an Algebraic Branching Program (ABP) computing a noncommu- 
tative polynomial in ¥{X) ( [N91) . |RS05j ). An ABP is a directed acyclic graph with one vertex 
of in-degree zero, which is called the source, and one vertex of out-degree zero, which is called the 
sink. The vertices of the graph are partitioned into levels numbered 0,1, ... ,d. Edges may only go 
from level i to level i + 1 for i = 0,1, . . . ,d—l. The source is the only vertex at level and the sink 
is the only vertex at level d. Each edge is labeled with a homogeneous linear form in the variables 
X. The size of the ABP is the number of vertices. 

The ABP computes a degree d homogeneous polynomial / G ¥{X) as follows. Fix any path 7 from 
source to sink with edges ei, 62, • • • , e^, where is the edge from level i — 1 to level i, and let £i 
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denote the linear form labelling edge Cj. We denote by the homogeneous degree d polynomial 
ii ■ £2 - ■ -^d (note that the order of multiplication is important). The polynomial / computed by 
the ABP is simply 

-yev 

where V is the set of all paths from the source to the sink. 

We will also consider a slight variant of the above definition where we allow multiple sources and 
sinks. A multi- output ABP P is defined exactly as above, except that we allow multiple sources at 
level and multiple sinks at level d. For each source s and sink t, an ABP Pgt may be obtained 
from P by removing all sources other than s and sinks other than t. Let S = {si, S2, . . . , Sq} and 
T = {ti, • • • ) ife} denote the sets of sources and sinks respectively in P. The ABP P will be 
thought of as computing the ab many polynomials computed by the ABPs Psi,tj- More precisely, 
the output of the ABP P is an a x 6 matrix A with entries from ¥{X) such that the (i, j)th entry 
of A, denoted A{i,j), is the polynomial computed by the ABP Ps^,tj- It is easily seen that we 
can write A as Ylm£Md{x) ^rn'm', where Am G F"^*; we will call A^ the coefficient matrix of the 
monomial m in the matrix A. 



3 The Hadamard Product 

A key notion we require for all our reductions is the Hadamard product of polynomials that was 
introduced in |AJSn9j . 

Definition 1. Given polynomials f,g(z ¥{X), their Hadamard product h = f o g is defined as 
follows: h is the unique polynomial in ¥{X) such that for any monomial m G JV[[X), the coefficient 
h{m) = f{m) ■ g{m). 

In [ AJS09t Theorem 5] we show that given a noncommutative circuit for polynomial / and an ABP 
for polynomial g we can efficiently compute a noncommutative circuit for their Hadamard product 
fog. However, the construction we present in |AJS09j modifies the noncommutative circuit for 
the polynomial /. Hence, it will not work if we are allowed only black-box access to /, which we 
require for certain applications in this paper. 

Suppose we have an efficient black-box algorithm for evaluating the polynomial / G ¥{X), where the 
variables in X take values in some matrix algebra (say, nxn matrices over a field F). Furthermore, 
suppose we have an explicit ABP for the polynomial g. Ideally, we would like to obtain an efficient 
algorithm for computing their Hadamard product fog over the same matrix algebra. 

However, what we can show is that we can put together the ABP and the black-box algorithm for 
/ to obtain an efficient algorithm that computes fog over F. This turns out to be sufficient to 
prove all our hardness results for the different noncommutative determinants. 

Theorem 2. Fix c? G N. Let Z = {zi, Z2, . . . , Zm} be a set of noncommuting variables and g G ¥{Z) 
be a homogeneous polynomial of degree d such that g is computed by an ABP P of size S. Then, 
there exist matrices Ai, A2, ■ ■ ■ , An G Ms{¥) such that for any homogeneous polynomial f G ¥{Z) 
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of degree d, fog = /(jdi^i, ^2-225 • • • j ^n-2n)(l, 5*). Moreover, given the ABP P, the matrices 
Ai, A2, . . . 1 An can he computed in time polynomial in the size of the description of P. 

Proof. Let the vertices of P be named 1, 2, . . . , S where 1 is the source of the ABP and 5 is the sink. 
Define the matrices Ai,A2, ■ ■ ■ , An G Ms(¥) as fohows: Ai{k, I) is the coefficient of the variable Zi 
in the hnear form labelhng the edge that goes from vertex k to vertex /; if there is no such edge, 
the entry Ai{k, I) = 0. For any monomial m = Zi^Zi^ ■ ■ ■ G A4d{Z), let A^ denote the matrix 
Ai^Ai^ ■ ■ ■ Ai^. We see that 

) ) (^«2 ) ■ ■ ■ i^id^id) 



Note that the coefficient g{m) of a monomial m = Zi^Zi^-'-Zi^ in g is just Ani{'i,S) = 
Z]fci,fc2,...,fcd-i6[S] 11^=1 where /cq = 1 and kd = S. Putting the above obser- 

vations together, we see that f{AiZi,A2Z2,...,AnZn){l,S) = Y^rneMd{Z) fi'^)'^rn{'i-, S)m = 
YlrneMdiz) f^^)9^^)^ — f ° 9- Since the entries of the matrices Ai, A2, ■ ■ ■ , An can be read 
off from the labels of P, it can be seen that Ai,A2,.-- ,An can be computed in polynomial time 
given the ABP P. This completes the proof. □ 

Remark 3. We note that the matrices Ai in the statement of Theorem\^ can actually he computed 
from the ABP even more efficiently, say, in uniform AC*^. 

The following corollary is immediate. 

Corollary 4. |AJS09j Given a noncommutative circuit of size S' for f E F(Z) and an ABP of size 
S for g G F(Z), we can efficiently compute a noncommutative circuit of size 0{S' S'^) for fog. 

The next corollary is the more useful version for this paper. 

Corollary 5. Let Z = {zi,Z2, ■ ■ ■ ,Zn}. Suppose A is a polynomial-time algorithm for computing 
a homogeneous degree d polynomial f G F(Z) for matrix inputs from Ms(F)Jl Given as input 
an ABP P, with S nodes, computing a homogeneous degree d polynomial g G F(Z), and scalars 
ai,a2, . . . , On G F, we can compute f o g[ai,a2, • • • , fln) in polynomial time. 

Proof. We first compute matrices Ai, A2, ■ ■ ■ , An, described in the Theorem [2l in time poly- 
nomial in the description of the ABP P. Then we invoke the given algorithm A on in- 
put {Aioi, A2a2, . . . , AnOn) to obtain as output an x 5 matrix whose (1,5')*'' entry contains 
/ o g{ai, 02, ... , On). Clearly, the simulation runs in polynomial time. □ 

^The statement can be generalized to any unital algebra A in place of the field F. 



f{AiZi,A2Z2,...,AnZn)= 22 fi^h^i2--- 

11,12,. ..,id&[n] 

~ / (^11^42 

11,12, •••>*dG[n] 

= ^ f{m)A„^m 
meMd(Z) 
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4 The hardness of the Cayley determinant 



We consider polynomials over an arbitrary field F (for the algorithmic results F is either rational 
numbers or a finite field). The main result of this section is that if there is a polynomial-time 
algorithm to compute the 2n x 2n Cayley determinant over inputs from Ms (F) for S = c - v? (for 
a suitable constant c) then there is a polynomial-time algorithm to compute the n x n permanent 
over F. 

Throughout this section let X denote {xij \ l <i,j < 2n}, and Y denote {yij | 1 < ^, j < n}. Our 
aim is to show that if there is a polynomial-time algorithm for computing Cdet2n(X) where Xij 
takes values in Ms{¥) then there is a polynomial-time algorithm that computes Cperm^(Y) where 
Uij takes values in F. 

The 2n x 2n determinant has 2n! many signed monomials of degree 2n of the form 
3^i,ct(i)^2,ct(2) ' " X2n,a{2n) a G S2n- We will identify n! of these monomials, all of which have 
the same sign. More precisely, we will design a small ABP with which we will be able to pick out 
these n! monomials of the same sign. 

We now define these n! many permutations from S2n which have the same sign and the correspond- 
ing monomials of Cdet2n that can be picked out by a small ABP. 

Definition 6. Let n G N. For each permutation tt G 5„, we define a permutation p{-k) in S2n, 
called the interleaving of tt, as follows: 



That is, the elements p(7r)(l), /3(7r)(2), • • • ,p{7r)(2n) are simply 7r(l),(n + 7r(l)), 7r(2), (n -|- 
7r(2)),--- ,7r(n),(n + 7r(n)). 

The following lemma states a crucial property of the permutation piir). 

Lemma 7. The sign of the permutation p(7r) is independent ofn. More precisely, for every vr G 5„, 
we have sgn(p(7r)) = sgn(p(l„)), where In denotes the identity permutation in Sn- 

Proof. For each vr G S'n we can define the permutation tt2 G S2n as 7r2(i) = 7r(i) for 1 < i < n 
and TT2{n + j) = n + 7r(j) for 1 < j < n. It is easy to verify that sgn(7r2) = sgn(7r)^ = 1 for every 
TT & Sn. To see this we write n2 as a product of disjoint cycles and notice that every cycle occurs 
an even number of times. Furthermore, we can check that p{Tr) = /9(l„)7r2, where we evaluate 
products of permutations from left to right. Hence it follows that sgn(p(7r)) = sgn(/9(l„))sgn(7r2) = 



We will denote by po the permutation p(ln), where 1„ denotes the identity permutation in S, 




sgn(p(l„)). 



□ 



For a G S2n, we will denote by nio- the monomial aJi,o-(i)^2,(7(2) " ' ' X2n,a{2n) ^ M{X). For a,T E S2n, 
we win denote the monomial X(^(i)^^(i)X(^(2),r(2) • • • Xa{2n),T(2n) by ma,T- 
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In the next lemma we show that there is an ABP that will filter out monomials that are not of the 
form ?7ip(7r) from among the m^. 

Lemma 8. There is an ABP P of size 0{in?) and width n that computes a homogeneous polynomial 
F E F(X) of degree 2n such that for any o", r G S2n, 

• F{mu) = 1 if a = p{tt) for some vr € Sn, and otherwise. 

• F{ma^r) = unless a = l2n, where l2n denotes the identity permutation in S2n- 
Moreover, the above ABP P can he computed in time poly(n). 

Proof. The ABP is essentially just a finite automaton over the alphabet X with the following 
properties: for input monomials of the form mo- it accepts only those monomials that are of the 
form mp^^-j . Further, for input monomials of the form mo-,r it accepts only those monomials of the 
form mi^^^r- We give the formal description of this ABP P below. 

The ABP P contains 2n + \ layers, labelled {0, 1, . . . , 2n}. For each even i S {0, 1, . . . , 2n}, there 
is exactly one node qi at level i; for each odd z S {0, 1, . . . , 2n}, there are n nodes Pi,i,Pi,2, ■ ■ ■ ,Pi,n 
at level i. We now describe the edges of P: for each even i £ {0, 1, . . . , 2n — 2} and j G [n], there 
is an edge from qi to Pi+i,j labelled Xj+ij; for each odd i G {0, 1, . . . , 2n} and j G [n], there is an 
edge from pij to gj+i labelled Xj+i^n+j. 

It is easy to see that P as defined above satisfies the requirements of the statement of the lemma. 
It is also clear that the ABP P can be computed in polynomial time. □ 

Note that the ABP P of Lemma [8] can in fact be constructed in uniform AC". 

Remark 9. For this section we require only the first part of Lemma\^ The second part of Lemma\^ 
is used in Section\^ 

We are now ready to prove that if there is a small noncommutative arithmetic circuit that computes 
the Cayley determinant polynomial, then there is a small noncommutative arithmetic circuit that 
computes the Cayley permanent polynomial. 

Theorem 10. For any n G N, if there is a circuit C of size s computing Cdet2n(^)) then there is 
a circuit C of size polynomial in s and n that computes Cperm„(y). 

Proof. Assuming the existence of the circuit C as stated above, by Corollary SI there is a noncom- 
mutative arithmetic circuit C" of size poly(s,n) that computes the polynomial F" = Cdet2n o F, 
where F is the polynomial referred to in Lemma El For any monomial m, if m 7^ mo- for any 
a G 5'2n, then Cdet2n('TT') = and hence, in this case, F"[m) = 0; moreover, for m = mo, we have 
F(m) = 1 if £7 = p(7r) for some tt G 5„, and otherwise. Hence, we see that 

= ^ sgn(p(7r))mp(^) = sgn(/3o) I ^ "ip(^) j 
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where the last equality follows from Lemma [71 



Let C be the circuit obtained from C" by substituting Xij with yi+j_ ■ if i is odd and j € [n], and 
by 1 if i is even or j ^ [n], and by multiplying the output of the resulting circuit by sgn(/9o). Let 
F' denote the polynomial computed by C . Then, we have 

p'i^) = E Km 

where "'-p(^) denotes the monomial obtained from mp(^^-^ after the substitution. It can be checked 
that for any vr € S'n, the monomial m'p(^) = 2/1,77(1)^2,77(2) '•• yn,7r(ra)- Hence, the polynomial F' 
computed by C in indeed Cperm„(y). It is easily seen that the size of C is poly(s,n). □ 

We now show that evaluating the polynomial Cdet2n over Ms{¥), for 5 = c • for suitable c> 0, 
is at least as hard as evaluating the permanent over F. 

Theorem 11. If there is a polynomial-time algorithm A that computes the 2n x 2n Cayley de- 
terminant of matrices with entries in Ms{F), for S = c ■ for suitable c > 0, then there is a 
polynomial-time algorithm that computes the n x n permanent over F. 

Proof. This is an easy consequence of Corollary [5l Consider the algorithm given by Corollary [5] 
for computing Cdet2n o F over the field F, where the ABP in Corollary [5] is the ABP of Lemma [8] 
computing F. 

In order to evaluate the permanent over inputs ajj, 1 < i, j < n we will substitute 2:24-1,^ = atj for 
^ ^ i, j ^ n- and we put Xij = 1 when i is even or j > n. As in the proof of TheoremfTUlit follows that 
for this substitution the algorithm computing Cdet2raoF will output sgn(yOo)Cperm„(aii, . . . ,ann)- 
Since sgn(po) can be easily computed, we have a polynomial-time algorithm for computing the nxn 
permanent over F. □ 

Remark 12. The above result has a stronger consequence: for any fixed e > 0, if there is a 
polynomial-time algorithm that computes the mxm Cayley determinant over Mms{F), then there 
is a polynomial-time algorithm that computes r2(m^/^) x r2(m^/^) permanents overY, hence implying 
that permanent over F is polynomial-time computable. 



5 The Cayley determinant over Clifford algebras 

We now consider the complexity of computing the determinant over real Clifford algebras of poly- 
nomially large dimension. We show via a polynomial-time reduction that computing the permanent 
over rationals is reducible to this problem. Indeed, by inspecting our result we can observe that 
even approximating the determinant over such Clifford algebras would yield similar approximation 
algorithms for the permanent over the reals. 

We first define the basic notions in the theory of Clifford algebras. A thorough treatment can be 
found in [LS09] . Fix m G N. The (real) Clifford algebra CL'^ is a 2™'-dimensional vector space 
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over M with basis elements of the form Ci^ei^ei^ ' ' ' ^ik where ii < 12 < is ■ ■ ■ < ik are elements 
from [m]. Multiplication between elements of the basis is defined by the following rules: = 1 
and CiCj = —ejCi for distinct i,j G [m]; this is extended linearly to all pairs of elements from 
the Clifford algebra. Given ii < 12 < ■■■ < ik from [m], we denote by 65 the basis element 
Cj^Cij ■ ■ ■ e-i^i where 5 = {ii, • • • > ^fc}- Each element of the Clifford algebra is uniquely expressible 
^ Ssc[m] ^565, where € M for each S. (Note that 60 and 1 both refer to the multiplicative 
identity of the algebra.) An idempotent of the Clifford algebra is an element e such that 6^ = e. 
Given h = J2sc[m] ^s^s in C-^m> "^^ define its norm \h\ to be ^^sci[m] "^l- 

The subset of basis elements {eg | S has even cardinality} generates a strict subalgebra of CL'„^. 
We will denote this subalgebra by CLm- This is the algebra of interest to us. The term 'Clifford 
algebra' will henceforth refer to CL^ for some m € N. 

Chien, Rasmussen, and Sinclair [CRS03] have shown that a polynomial-time algorithm that, when 
given as input an n x n matrix B with entries from CL^ for m = 21ogn + 2, computes |Cdet„(-B)|^ 
can be used to design a randomized polynomial time algorithm to approximate the 0-1 permanent 
(over Q). 

In this section, we prove that if there is a polynomial-time algorithm to compute either |Cdetri(-B)P 
or Cdet„(i?), then the permanent (over inputs from M) can actually be computed in polynomial 
time. For an n x n real matrix A, let perm„(A) denote the permanent of A. 

Remark 13. In a sense, our result in this section should not he surprising. We have already 
proved (in Theorem that computing the determinant over matrix algebras is at least as hard 
as computing the permanent. Also, it is known that Clifford algebras of polynomial dimension are 
isomorphic to matrix algebras of polynomial dimension (see, for example, JLSOS\ Chapter 5]). How- 
ever, in this section we actually give an explicit polynomial-time reduction showing that computing 
the permanent over the reals is reducible to computing either |Cdet„(S)p or Cdetn{B) where the 
entries of B are from the Clifford algebra CLm- 

Suppose we wish to compute the permanent of an n x n matrix with entries from R. W.l.o.g., we 
assume that n = 2^ for some £ G N. Let m denote 5£. The next lemma is about the existence of 
certain elements in the algebra CL^ useful for the reduction. 

Lemma 14. Let n,£,m be as above. Then, there exist hi,h2, ■ ■ ■ ,hn,h'ih2, ... ,h'^ G CLm and an 
idempotent e £ CLm such that: 

• For all j, hjh'j = e. 

• For all j 7^ k, hjh'^ = 0. 

I [2 1 

• |e| = 27- 

Moreover, the elements hi, h2, . . . , hn, h'l, h2, . . . , h'^ and e can be constructed in time poly(n). 
We defer the proof of the above lemma and first prove the main result of this section. 
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Theorem 15. Let n,£,m be as above. There is a polynomial-time algorithm which, when given 
any matrix A E M„(M), computes a B £ M2n{CLm) such that |Cdet2n(5)|^ = ^^^"^^/-^^ . 

Proof. The matrix B wih be the following: for any odd i € [2n] and any j € [2n], set B{i,j) - the 
(i, j)th entry of - to be A{^^,j)hj if j < n and if j > n; for any even i G [2n] and j S [2n], 
set B{i,j) to be /ij_„ if j > n and otherwise. Clearly, B can be computed in polynomial time 
given A. Note the following property of B: for any odd i G [2n] and j, k G [2n] 

Bii i)B(i + l k)-i ^(^'-^^ ifi<nand/>: = n + j, 
B[t,j)B[2 + i,k)-<^ ^ otherwise. 

Here e denotes the idempotent from Lemma O The following claim is easy to see. 

Claim 16. For any permutation a € S2n, the product Yli^i B{i,a{i)) = {YVi=i ^ih'^i'i')))^ if 
a = p{tt) for some tt £ Sn and it is otherwise. 

Let us consider Cdet2n(-B). We have: 

Cdet2n(B) = Yl sgn((7)5(l, a(l)) • B(2, a{2)) ■ ■ ■ B{2n, a{2n)) 

0"GS'2n 

n 

= J2 sgn{p{7r))([lA{i,TT{i)))e 

TT&Sn j = l 

= sgn(po)perm„(A)e 

Thus, we see that |Cdet2n(5)P = perm„(A)2 |ep = □ 
We have the following easy consequence of the above theorem. 

Corollary 17. Fix any e > 0, and suppose there is a polynomial-time algorithm that computes 
|Cdet„(-B)|^ on input an n x n matrix B with entries from CLm for m = elogn. Then there is a 
polynomial-time algorithm that computes the nxn permanent of matrices with nonnegative rational 
entries. 

Proof. The statement directly follows from Theorem [15] for m = [5 log n] . To prove hardness for 
m = elogn, we note that a polynomial-time algorithm to compute |Cdet„(i?)|^ over CL^iogn can 
be used to compute |Cdet„e/5(5)|^ over CLgj^g^e/s in polynomial time. □ 

A (5-approximation algorithm A for a function / : S* — > Q is an algorithm such that for each 

{l-6)f{x)<A{x) < (1 + 5) fix). 

Our reduction from computing the permanent for nonnegative entries to computing |Cdet„(-B)p 
actually yields an approximation preserving reduction. We formalize this in the next corollary. 
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Corollary 18. Fix any 6 > and e > 0. Suppose there is a polynomial-time 5 -approximation 
algorithm for the function that on input an n x n matrix B with entries from CLm for m = elogn 
takes the value |Cdetn(-B)p. Then there is a polynomial-time 5 -approximation algorithm for the 
n X n permanent with nonnegative rational entries. 



We now prove Lemma [T^ 



Proof of Lemma [74 Let ei, 62, . . . , denote the generators of CL^. Partition the set [m] into £ 



subsets of size 5 as foUows: set Si = {5{i — 1) + i | j G [5]} for each z G For each i G [^], let 
^.^o = {5(i - 1) + 1, 5(i - 1) + 2, 5{i - 1) + 3, 5{i - 1) + 5} and 5^,1 = {5{i - 1) + 2, 5{i - 1) + 3, 5{i - 
l) + 4,5(i-l) + 5}. 

Using the fact that e? = 1 and ejCj = —ejei for i ^ j it easily follows that for any two disjoint 
sets S, r C [m] such that \S\, \T\ are even, we have eser = exes- Hence, the elements 65-. and 
eSj commute for i ^ j and any 61, 62 G {0, 1}. Furthermore, for all i S [i] and b € {0, 1} we have 
e'g, ^ = 1. Also, we have 65. qCs- ^ = —65. ^e^- q. Finally, notice that 65^ for 1 < i < £ and b S {0, 1} 
are all elements of CLm- 



For i € [£] and b G {0,1}, set gifi = — ^'^d gi^i = — — Also, set g'iQ = gi.o and 

g'ii = 2 Notice that gf^ = gi^. We also note an additional relation 65. (,(1 — £5- j) = 

(1 + 65. 1)65. Q. Using these we can easily derive the following crucial properties of these elements 
of CLm- 

• For each i £ [i] and 6 € {0, 1}, = difl- 

• For each i £ [i] and b G {0, 1}, gi^b9i i-b ~ 

• For ii ^ i2 and any 61, 62 £ {0, 1}, the elements gi^fii and g-^ commute. 

Finally, we define hj, h'j for a fixed j S [n]. Let 5i 62 ... 6^ be the binary representation of the integer 
j - I (recall that n = 2^). We define hj = 51,61 92,62 ' ' ' 9e,be and /i^- = g'i^hi92,b2 ' " dtbg- Also, we 
define e to be 91,092,0 • • • 9^,0) which is the same as hi and h[. 

We now prove that the hj, h'j {j S [n]) and e satisfy the properties claimed in the statement of the 
lemma. Fix any j G [n] and let bib2 - ■ - bi be the binary representation of j — 1. We have 

hjh'j = 9i,fei 92,62 • ■ ■ 9i,be9i,bi92,b2 ■ ■ ■ 9^,6^ 

= (91,61 9l,6i) • (92,6192,62) ■ • • {9i,be9l,be) 
= 91,092,0 ■ • • 9^,0 = e 

The second equality follows from the fact that gi^^^b and g'- commute for any distinct ii and 12- 
The third equality follows from the fact that for any i and b, gifig[ ^ = gi^- This proves the first 
property claimed in the statement of the lemma. Similarly, we can see that e is an idempotent: 
e^ = h\ = e. 
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Fix any distinct j, k G [n]. Let 6162 ■ ■ - be and 6'^ 63 ■ ■ - b'i be the binary representations of j and k. 
Since j 7^ A:, we can fix some i such that 6j ^ 6^. We have 

hjK = 9iM92,b2 ■ ■ ■ 9e,bt9i,b[92,b'2 ' " " 9^^^ 

= i9iM9'iM^ ■ (52,6152,62) • • • (.9iM9i,b'^ ■ ' ' (.9e,be9e,bJ 

= {9iM9i,bJ ■ (52,61 52,62) • • • • • • {gi,be9e,bt) = 

where the third equahty follows from the fact that we have gi^b9i i_b = 0. This proves the second 
claim made in the lemma. 



Finally, we note that 



|ep = j5i,o52,o • • •5£,oP 



1 



TceieT 



1 



It is easily seen from their definitions that the hj, h'j and e can be computed in time poly(n). This 
completes the proof of the lemma. □ 



6 The Symmetrized Determinant 

In this section, we observe that the 2n x 2n symmetrized determinant over 0(n^)-dimensional matrix 
algebras is at least as hard to compute as the permanent. This stands in marked contrast to the 
result of Barvinok [B], who shows that over constant- dimensional matrix algebras, the symmetrized 
determinant is polynomial-time computable. 

In this section, let F denote a field of characteristic 0. Let X = {xij | 1 < i,j < 2n} and Y = 
{Vij I 1 < i,j < n}. Recall that for cj, r G 52„, the monomial m^,^ is x^(i)^^(i)X^(2),r(2) ■ ■ ■ Xa(2n),T{2n), 
and the monomial m„ is Xi^^(i)X2,<j(2) • • • X2n,a{2n)- 

Theorem 19. If the sdet2n(^) polynomial over¥ can be computed by a polynomial- sized noncom- 
mutative arithmetic circuit, then the polynomial Cperm„(y) can also be computed by a polynomial- 
sized noncommutative arithmetic circuit. 

Proof. Assume that sdet2n(-'^) is computed by a circuit C of size s. As in Theorem [TOl we will 
proceed by taking Hadamard product. Let P be the ABP defined in Lemma [8] and F(X) the 
polynomial it computes. Let F" denote the polynomial sdet2n(-^) ° F. Note that by Corollary U 
F" can be computed by a circuit C" of size poly(,s, n). From LemmaEl we have FlrUa^r) = unless 
o" = l2n) the identity permutation in S2n', moreover, we also have F[mi,^^^r) = F{mr) which is 1 if 
r = p('/r) for some n £ Sn and otherwise. By the above reasoning, 

^"^^^ " (i)! E sgn(/>(7r))mp(,) = ^-^^ E "^p(-) 
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Now, we substitute each Xij by yi+i ■ if i is odd and j E [n] and by 1 if z is even or j ^ [n] in 
the circuit C". The effect of this substitution is to transform mpj-j^) into yi,7r(i)2/2,7r(2) ' ' ' Vn,-K{n) for 
each vr € S'n- Hence, the resulting polynomial is simply '^g'^(Po)Cpcrm„(y) ^ ^pj^^g^ )-,y multiplying by 
sgn(po)(2n)!, we obtain a circuit C of size poly(s,n) that computes Cperm„(y). □ 

Theorem 20. If there is a polynomial-time algorithm A that computes the 2n x 2n symmetrized 
determinant of matrices with entries in Ms{¥), for S = c ■ for suitable c > 0, then there is a 
polynomial-time algorithm that computes the n x n permanent over ¥. 



Proof. The proof is almost exactly identical to that of Theorem [TTJ Consider the algorithm given 
by Corollary [5] for computing sdet2n o F over the field F, where the ABP in Corollary [5] is the ABP 
of Lemma [8] computing F. 

In order to evaluate the permanent over inputs ajj, 1 < i,j < n we will substitute X2i—\^j — 
for 1 < i,j < n and we put Xij = 1 when i is even or j > n. As in the proof of The- 
orem [T9l it follows that for this substitution the algorithm computing sdet2n o F will output 
^^^n)\^ Cperm^(aii, . . . ,a„„). Since sgn(po) and (2n)! are easily computable, we have a polynomial- 
time algorithm for computing the n x n permanent over F. □ 



7 The Moore determinant 



We demonstrate by a simple reduction that the Moore determinant and permanent are interre- 
ducible. We also show that the computing the Moore determinant over a field of characteristic 
zero is at least as hard as counting the number of directed Hamilton Cycles of a directed graph, 
which is a well-known ^^P-complete problem. If the field is of characteristic k, then computing the 
Moore determinant over the field is at least as hard as counting the number of Hamilton cycles of 
a directed graph modulo the prime k, which is hard for Mod^P. 

Assume X = {xij | 1 < i,i < n}. Given a permutation a € Sn, we write cr as a product of 
disjoint cycles as follows: {nf^ ■ ■ ■ f^i;J(f^2i ' " " ''^212^ ' ' ' ("-ri ' ' ' "-w^) '^^^^ "-fi < ^ij ^ ^ V] 

and j G \ {1} and satisfying > rigi > •■■ > n'^i- Let denote the monomial 

Xri'^ ' ' ' Xri^ ' ' ' Xri^ ' ' ' Xri^ • 

nil,ni2 '^hlr' 11 "'rl'"'r2 "'rlr' r\ 

Let Cn denote the set of all 1-cycles in Sn, i-e permutations whose cycle decomposition con- 
sists of a single cycle of length n. Define the polynomial HC„(xii, . . . , G ¥{X) to be 
So-eC„ ^o"- ^'^y directed graph G on n vertices with adjacency matrix A. Let H[G) denote 
HC„(^(1, 1), . . . , A(n, n)). The quantity H{G) has a simple description: if F is of characteristic 0, 
then H{G) is the number of directed Hamiltonian cycles in G; and if F is of characteristic k, then 
H{G) is the number of directed Hamiltonian cycles of G modulo k. 

We have the following easy lemma: 

Lemma 21. There are ABPs P[ and P2 of size 0{n'^) and width n that compute homogeneous 
polynomials F{,F2 G ¥{X) of degree n such that for any a G Sn, we have 
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• F{{w^) = sgn(o-). 



• F2{wa) = sgn((T) if a G Cn and otherwise. 
Moreover, the above ABPs can be computed in time poly(n). 

Proof. Recall that given a permutation a G 5„, the quantity sgn(fT) is (—1)""'"'^'^, where is 
the number of disjoint cycles in a. Moreover, note that if cr as a product of disjoint cycles is 
(nJi • • • '^iij)('^2i ■ ■ ■ ''^212^ ' ' ' (""ri ■ ■ ■ "-w^) ^ above, the value is simply the number of left-to-right 
minima in this representation, i.e the number of n^. such that n^. < n^^ for all n^^ to the left of 
nfj. Using this observation, it is easy to design an ABP P[ that keeps track of the sign of the 
permutation and computes a polynomial F[ as above. The ABP can be constructed similarly; 
the main difference from the case of P{ is that the ABP must produce the coefficient unless 
n'^i = 1. We omit the formal description of P[ and I— ' 

The analogue of Theorem [10] for the Moore determinant follows below. The statement here is 
stronger: we show that the arithmetic circuit complexity of Mdet„(X) is polynomial if and only if 
the arithmetic circuit complexity of Mperm„(X) is polynomial. 

Theorem 22. The Moore determinant polynomial Mdet„(X) can be computed by a polynomial- 
sized noncommutative arithmetic circuit if and only if the Moore permanent polynomial Mperm„(X) 
can be computed by a polynomial- sized noncommutative arithmetic circuit. 

Proof. As in the proof of Theorem 1101 we will use the Hadamard product; this time, it can be used 
to erase or introduce the signs of the permutations corresponding to each monomial w^. Formally, 
we have Mperm„(X) = Mdet„(X) o F[{X) and Mdet„(X) = Mperm„(X) o F[{X), where F[{X) is 
the polynomial defined in the statement of Lemma[2TJ Hence, if Mdet„(X) (resp. Mperm„(X)) is 
computed by a noncommutative arithmetic circuit of size s, then by applying Corollary [U we see 
that Mperm„(X) (resp. Mdet„(Ar)) is computed by a noncommutative arithmetic circuit of size 
poly(s,n). □ 

Remark 23. Note that Theorem\2^ proves an equivalence (up to polynomial factors) between the 
arithmetic circuit complexities of the Moore determinant and permanent. This is a stronger state- 
ment than we obtained in the case of the Cayley determinant and permanent, where we only showed 
(roughly) that the Cayley determinant is at least as hard to compute as the Cayley permanent. The 
reason for this is that we are unable to obtain a small ABP that performs the function of P[ for the 
monomials m„ (defined in Section\^: that is, a small ABP computing a polynomial Fi such that 
Fi^niu) = sgn(mo-) for every a £ Sn. However, we are unable to rule out the possibility that such an 
ABP exists. If it does, then as above, we can obtain a simple equivalence between the complexities 
of the Cayley determinant and permanent. 

We now consider the complexity of computing the Moore determinant over matrix algebras of 
polynomial dimension. We can, as in the previous sections, show that this is at least as hard as 
computing the permanent over matrices with entries from F, but we take a different route this 
time. We show that if the Moore determinant over a field of characteristic k can be computed 
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in polynomial time, then there is a polynomial-time algorithm to compute the number of directed 
Hamilton cycles H{G) modulo k for an input directed graph G. This allows us to draw stronger 
consequences, namely that the Moore determinant is hard to compute even when the field F is 
of characteristic 2, something that would not follow if we reduced the permanent to this problem 
(since the permanent is polynomial-time computable over fields of characteristic 2). 

Theorem 24. If there is a polynomial-time algorithm A that computes the nxn Moore determinant 
of matrices with entries in ^^^(F), for S = c-n^ for suitable c > 0, then there is a polynomial-time 
algorithm that, on input a directed graph G, computes H{G). 

Proof. Note that HC„(X) = Mdet„(X) o F^, where is the polynomial computed by ABP 
constructed in Lemma [211 Moreover, H{G) = HCn(A(l, 1), . . . , A{n, n)), where A is the adjacency 
matrix of the graph G. Hence, to compute H{G), we need to compute HC„(A(1, 1), . . . ,A{n,n)), 
which can be done in polynomial time by Corollary [5l □ 



8 Completeness Results 

In this section we observe that the noncommutative Cayley determinant over integer matrices 
is complete for GapP w.r.t. polynomial-time Turing reductions. Likewise, the noncommutative 
Cayley determinant over a finite field of characteristic A; 7^ 2 is hard for the modular counting 
complexity class Mod^P w.r.t. polynomial-time Turing reductions. These observations also hold 
for the symmetrized determinant. For the Moore determinant, we prove the above results without 
any restriction on the characteristic of the underlying field. We formally describe these observations. 

Definition 25. [FFK94j . [BG92j A function f : S* — > Z is in GapP if there is a polynomial time 
NDTM M such that for each x G S* the value f{x) is a.ccM{x) — Te]i^/[{x). 

For a prime k, the class Mod^P consist of languages L C S* such that for some function f G GapP 
we have x (z L if and only if f{x) = 0{mod k). 

By Valiant's result [V79j it is known that the integer permanent is GapP-complete with respect to 
polynomial-time Turing reductions. Furthermore, the permanent over F^ is Mod^P-hard for prime 

Now, for n G N, consider the Cayley determinant for 2n x 2n matrices with entries from MsC^), 
where S = cn^ for some constant c. By Theorem [TTl there is a fixed c > such that computing the 
integer permanent for nxn matrices is polynomial-time reducible to computing the (1, SY^ entry 
of such a Cayley determinant. The same observation holds modulo k for a prime k. 

Furthermore, the problem of computing the (1, 5)*^ entry of such a Cayley determinant over Z is 
easily seen to be in GapP: we can design a polynomial-time NDTM which takes as input a 2n x 2n 
matrix with entries from M5(Z) and the difference in the number of accepting and rejecting paths 
is the (1, entry of its Cayley determinant. Hence we have the following. 
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Corollary 26. There exists a constant c such that the following holds. For S = cn^ , computing the 
(1,5)*^ entry of the Cayley determinant for 2n x 2n matrices with entries from Msi"^) is GapP- 
complete w.r.t. polynomial-time Turing reductions. Given a finite field F of characteristic k ^ 2, 
computing the (1,5)*'' of the Cayley determinant for 2n x 2n matrices over Ms{¥) is hard w.r.t. 
polynomial-time Turing reductions for ModfcP. 

We have similar GapP-completeness and Mod^P-hardness consequences for the symmetrized deter- 
minant from the results in Sections [H For the Moore determinant, by Theorem 1241 we additionally 
obtain hardness for 0P over fields of characteristic 2. 

Corollary 27. There exists a constant c such that the following holds. For S = cn^ , computing 
the (1, 5)*^* entry of the Moore determinant for 2n x 2n matrices with entries from MsC^) is GapP- 
complete w.r.t. polynomial-time Turing reductions. Given a finite field F of any characteristic 
k > 1, computing the {1,SY^ of the Moore determinant for 2n x 2n matrices over Ms(F) is hard 
w.r.t. polynomial-time Turing reductions for ModkP . 

Proof. The result follows from Theorem 1241 and the following observations: computing H{G) over 
the rationals on an input graph G is GapP-complete w.r.t. polynomial-time Turing reductions; 
similarly, computing H[G) over a field F of characteristic k (including k = 2) \s hard for Mod^P 
w.r.t. polynomial-time Turing reductions. □ 

9 Discussion 

Our work raises further interesting questions regarding the complexity of the noncommutative 
determinant. 

An important open question is the complexity of computing the noncommutative determinant over 
constant dimensional matrix algebras. Theorem [11] can be easily used to show that assuming that 
the permanent of an n x n matrix over F cannot be computed in subexponential time, the n x n 
noncommutative Cayley, symmetrized, and Moore determinants with entries from M^2ogn)"(i) (^) 
cannot be computed in polynomial time. Can one strengthen this result to one that says something 
about computing the Cayley or Moore determinant over matrices with entries from Mc(F) for some 
absolute constant c? (Recall that the symmetrized determinant, on the other hand, is efficiently 
computable over constant dimensional matrix algebras.) It is interesting to note that |CS04| have 
shown an exponential lower bound for the ABP complexity of the Cayley determinant over even 
2x2 matrices. 

A question that arises from the results of Section H] is if one can show - analogous to those results 
- that the Cayley permanent is at least as hard to compute as the Cayley determinant. As pointed 
out in Remark 1231 one way to prove this is to construct a small ABP Pi that computes a polynomial 
Fi of degree n such that ^(2:1 o-(i) • • • Xn^a{n)) = sgn{a) for every a G Sn- This would also make the 
proofs of Theorems 1101 and 1111 much more transparent. 

Finally, note that our results do not imply that the Cayley determinant is hard to compute over 
Mfc(F) when F is a field of characteristic 2, since the permanent is known to be polynomial-time 
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computable over such fields. On the other hand, we have proved that the Moore determinant over 
such domains (where k is polynomial) is hard for ©P. Can we prove an analogous result for the 
Cayley determinant? 
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